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Numerical solutions of the Wong- Yang-Mills equations with anisotropic particle mo- 
mentum distributions are presented. Their isotropization by collective effects due to the 
\ classical Yang-Mills field is shown. 
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1. INTRODUCTION 

It is of great importance to understand the non-equilibrium dynamics in the early stage 
of high energy heavy-ion collisions and if, how, and when the produced partons thermalize. 
Boltzmann transport approaches with inelastic processes have been employed to follow 
the scattering of the hard gluons among themselves [QJ. However, the soft classical color 
+f | fields may give rise to non-perturbative collective processes and hence could also play an 
'"q ■ important role. Specifically, QCD plasma instabilities may develop, e.g. due to anisotropic 
distributions of released hard partons As a consequence, the "bottom-up scenario" [ 
|3| for thermalization of the Color Glass Condensate would be modified significantly [ H] . 
.£h [ Indeed, within the hard loop effective theory unstable soft modes were found if the 
hard modes exhibit an anisotropic distribution in momentum space. The time evolution 
has been studied via the full non-linear hard loop effective action [Ej. 3 + ID simulations 
of that theory [|7j show that instabilities grow more slowly once the field strength becomes 
large and non-Abelian self-interactions become nonperturbative. Possible effective poten- 
tials for anisotropic QCD plasmas beyond the hard loop approximation are discussed in 
Ref. [E|. 

In the present work we follow numerically the evolution of the hard modes, represented 
by particles, coupled to a classical SU(2) Yang-Mills field. The back-reaction of the 
fields on the particles is taken into account fully. We present some additional results as 
compared to ref. such as the modification of the particle spectra. 

2. Model 

We solve Wong's equations [ ITU] 

^ = » <) ^ = gQUE a + v t x B a ) , ^^igv^Ql (1) 
for the i-th (test) particle, coupled to the Yang-Mills equation 

D fl F^=j" = g y £Q^'6(x-x i ). (2) 
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These equations reproduce the "hard thermal loop" effective action [El near equilibrium. 
Numerical techniques to solve the classical field equations coupled to particles have been 
developed in Ref. [H2]- Our update algorithm is closely related to the one explained there. 

In the following, we assume that the fields only depend on time and on one spatial 
coordinate, x, which reduces the Yang-Mills equations to 1+1 dimensions. The parti- 
cles are allowed to propagate in three spatial dimensions. For simplicity, we also re- 
strict ourselves to the case without expansion here; the more realistic case with lon- 
gitudinal expansion (solved without particles in [ IT3]) will be addressed in the future. 
The initial anisotropic phase-space distribution of hard gluons is taken to be f(p,x) oc 
exp(—^Jpy + p^/phard) 5(p x ) . This represents a quasi-thermal distribution in two dimen- 
sions, with "temperature" = Phard- It is assumed to come about by free streaming of the 
particles for a time of order ~ 1/phard after the collision of the nuclei. In a comoving 
frame (with the central rapidity region) then, the longitudinal momenta are red-shifted 
to much below the typical transverse momentum. 

The initial field amplitudes are sampled from a Gaussian distribution with a width 
tuned to a given initial energy density. We solve the Yang-Mills equations in A = 
gauge and also set A = (i.e. all gauge links =1) at time t = 0; the initial electric field 
is taken to be polarized in a random direction transverse to the x-axis. 

3. RESULTS 
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Figure 1. Time evolution of the pressures for Figure 2. Time evolution of the pressures 
weak initial fields. for the strong field case. 



Fig.^shows the time evolution of the particle and field pressures for the initial condition 
Phard = 10 GeV and initial field energy density ~ 10~ 2 GeV/fm 3 /<? 2 , for which instabilities 
have been seen in Ref. [ E] . We observe a rapid exponential growth of the longitudinal 
pressure, approaching an isotropic configuration. This is due to deflection of the particles 
in the exponentially growing field, which randomizes their momenta. Note also that the 
dominant contribution to the build-up of longitudinal pressure is from deflection of the 
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particles, not from the fields themselves. In our 1+1D simulations the transverse pressure 
of the fields is zero T^ ld = 0. Furthermore, we find that T^ ld = Tg£ d , i.e. equal energy 
density and pressure, which is due to the dominance of transverse magnetic fields. Since 
periodic boundary conditions are used, no collective flow is generated: T° l = 0. 

The separation between hard and soft modes may not be very large at RHIC energy [ 
IT3] . In Fig |2l we show the time evolution of the pressures for Phard = 1 GeV and initial 
field energy density ~ 10 _1 GeV/fm 3 /g 2 ; the time scale is set by the lattice size L, =10 fm 
for this simulation. When the separation between hard and soft modes is small, strong 
instabilities are not seen, but the system still approaches isotropy [E]. In fact, the rapid 
isotropization of the particle momenta by the strong fields prevents the occurrence of 
pronounced instabilities. 
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Figure 3. Particle distribution functions for 
the weak field case. 



Figure 4. Particle distribution functions 
for the strong field case. 



It is interesting to look at the particle distribution function itself, even though we 
presently restrict ourselves to simulations of collisionless plasmas. Field fluctuations gen- 
erate an effective collision term [E] mediating soft exchanges. We plot initial as well as 
final (t/L/N c = 0.5) particle distribution functions for weak initial field in Fig. |3J and 
for strong field in Fig. |3J For the former case, we observe that the initial ^-function in 
the longitudinal direction becomes a thermal distribution extending over several orders of 
magnitude. However, the longitudinal direction is still much colder than the transverse 
directions, in agreement with the behavior of the pressures above. This is because of 
the assumption of transversally homogeneous fields. The fluctuations of the gauge field 
induce fluctuations in the motion of the quasiparticles. The presence of fluctuations is 
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closely linked to dissipative processes. 

In the future, full 3D simulations of the Wong- Yang-Mills equations are required. They 
will also allow us to study the effects due to a (azimuthally asymmetric) boundary. We 
also need to include longitudinal expansion of the system to check whether the distribution 
stays isotropic. Finally, collisions with momentum transfer above that mediated by the 
classical field have to be included for a complete description of the non-perturbative 
plasma in weak coupling. 
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